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Abstract 

The Casimir force due to thermal fluctuations (or pseudo-Casimir force) 
was previously calculated for the perfect Bose gas in the slab geom- 
etry for various boundary conditions. The Casimir pressure due to 
quantum fluctuations in a weakly-interacting dilute Bose-Einstein con- 
densate (BEC) confined to a parallel plate geometry was recently cal- 
culated for Dirichlct boundary conditions. In this paper we calculate 
the Casimir energy and pressure due to quantum fluctuations in a zero- 
temperature homogeneous weakly-interacting dilute BEC confined to a 
parallel plate geometry with periodic boundary conditions and include 
higher-order corrections which we refer to as Bogoliubov corrections. 
The leading order term is identified as the Casimir energy of a mass- 
less scalar field moving with wave velocity equal to the speed of sound 
in the BEC. We then obtain the leading order Casimir pressure in a 
general three-dimensional rectangular cavity of arbitrary lengths and 
obtain the finite-size correction to the parallel plate scenario. 
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1 Introduction 



After nearly 50 years since its prediction in 1948 by Casimir PP, the Casimir 
force has now been successfully measured by a modern series of experiments 
starting with Lamoreaux's 1997 landmark experiment ^ with a torsion pen- 
dulum which reduced errors dramatically compared to the early 1958 exper- 
iment by Spaarnay [Hj • The force was subsequently measured more precisely 
in 1998 using an atomic force microscope 11 and the measurements agreed 
with theoretical predictions to within 1% after finite conductivity, roughness 
and temperature corrections were taken into account. The Casimir force in 
the more difficult parallel plate geometry was then successfully measured to 
15% accuracy 5_. Thus the modern era of precise Casimir measurements was 
born. All the measurements of the Casimir force to date have been limited to 
the case of the electromagnetic field. However, experiments may soon mea- 
sure the Casimir force for a massless scalar field via quantum fluctuations 
in Bose-Einstein condensates (BEC). It has been noted (see ^;6J that the 
quasiparticle vacuum in a zero-temperature dilute weakly-interacting BEC 
should give rise to a measurable Casimir force. The fact that a scalar field 
propagates at the speed of sound in a BEC medium in contrast to the speed 
of light in Minkowski spacetime does not change anything fundamental in 
relation to the Casimir energy. If the speed of propagation is constant in a 
given medium, the Casimir energy in units of this speed will be the same 
value regardless of whether the medium is spacetime or a BEC. Moreover, 
a generally covariant action analogous to what we see in General Relativity 
exists for scalar fields propagating in a particular fluid. The Lagrangian is 
similar to that of a massless Klein-Gordon field with the Minkowski metric 
r]^ y of spacetime replaced by an effective or acoustic metric y 8, • Quoting 
directly from [S] , "at low momenta linearized excitations of the phase of the 
condensate wavefunction obey a (3-|-l)-dimensional d'Alembertian equation 
coupling to a (3-|-l)-dimensional Lorentzian-signature 'effective metric' that 
is generic, and depends algebraically on the background field.". 

The Casimir force due to thermal fluctuations in a perfect Bose gas confined 
to a slab geometry was recently calculated for Dirichlet, Neumann and pe- 
riodic boundary conditions 1(7 (see also comment on the work of [117 ) . 
The Casimir pressure due to quantum fluctuations for a weakly-interacting 
dilute Bose-Einstein condensate confined to a parallel plate geometry was 
also recently calculated [Jj. In this paper, we extend the work of [Jj on 
Casimir forces in BEC's to include a general three-dimensional cavity and 
consider periodic instead of Dirichlet boundary conditions. We first calcu- 
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late the Casimir pressure due to quantum fluctuations of the quasiparticle 
vacuum in a zero-temperature homogeneous dilute weakly-interacting BEC 
confined to a "parallel plate" geometry with periodic instead of Dirichlet 
boundary conditions. We show that the leading order term for the Casimir 
energy is equal to that of a massless scalar field moving with wave velocity 
equal to the speed of sound in the BEC. We also obtain the much smaller 
'Bogoliubov' corrections due to the nonlinearity of the Bogoliubov dispersion 
relation. We then generalize the results to a three-dimensional "rectangular" 
cavity of arbitrary lengths and obtain the leading order finite size correc- 
tions to the parallel plate scenario for periodic boundary conditions. The 
quotes around "rectangular" or "parallel plate" are simply to remind the 
reader that with periodic boundary conditions we are of course dealing with 
a hypertoroidal geometry. We drop the quotes from now on. 

To maintain a homogeneous gas in a slab geometry requires periodic bound- 
ary conditions. Though such boundary conditions constitute an idealiza- 
tion, it enables one to obtain relevant analytical results for the Casimir 
pressure of the BEC prior to a numerical analysis of the Casimir force in 
a non-homogeneous gas confined to a disk-like geometry via an anisotropic 
harmonic potential. 



2 Casimir pressure for BEC in parallel plate ge- 
ometry and massless scalar fields 

Consider a weakly interacting BEC characterized by an interparticle contact 
pseudopotential, 8 vr a (5^^^ (r), where a is the 2-particle positive scattering 
length (we work in units of 7i = 1 and 2m = 1). For this delta function 
potential, the contribution E of the depletion to the ground state energy 
due to quantum fluctuations in a zero-temperature untrapped homogeneous 
dilute weakly-interacting BEC is given by |12| ll.Sj 

k^O k^O 

where /u = 8 vr a p, is the number of atoms, p is the density N/V and 

E{k) = Eb - k'^ - = k \/A;2 + 2fi-k^ - p.. (2) 
Eb = k ^/k'^ + 2 n is the Bogoliubov dispersion relation. 
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With periodic boundary conditions, the homogeneous gas is confined to a 
parahel plate geometry with a trapping potential which is zero everywhere. 
It is worth noting that for quantum fluctuations to manifest themselves, 
the plate separation d must be much greater than the healing length i.e. 
d» /i^-"^/^. In an Li X L2 X d hypertoroidal space, k'^ = (2n7r/d)^ + 
(2ni7r/Li)^ + (2n2 7r/L2)^ where (n,ni,n2) 7^ (0,0,0) and n, ni and n2 
are integers that run from —00 to +00. Both the volume F = Li L2 d and 
the number of atoms are assumed large with the density p = N/V low 
enough that the dilute condition \J pa?" « 1 is satisfied. The triple sum for 
the ground state energy can be broken up in the following convenient 
fashion: 

00 00 00 

^Eik) = E E E ^(^)+^ 

k^O ra=— 00 ni=— 00 n2=— 00 

00 00 00 00 00 

= 8^^^ ^S(fc) + 8 5^ 5]i^(fc)+8/i. (3) 

n=l 711=0 n2=0 ni=0n2=0 



For parallel plates, Li and L2 are very large (infinite limit) and the sums over 
rii and 722 become integrals. The double sum in @ becomes a double integral 
and does not contribute to the Casimir energy being purely a continuous 
term. The constant 8^ also does not contribute to the Casimir energy. The 
relevant term for the Casimir energy, the triple sum in (jH)), becomes a sum 
over a double integral and k'^ reduces to (2 n vr/d)^ + r^. The relevant energy 
-E per unit area is then given by 




(5) 

The integral can be evaluated but this is not necessary since only the deriva- 
tives of f{n) will be needed to determine the Casimir energy. The parameter 
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A is an ultraviolet cut-off introduced because the energy E given by and 
hence f{n) are formally divergent. There is nothing physical about this di- 
vergence. It simply reflects that the delta function contact pseudopotential 
cannot be naively extrapolated to very high momentum where the wave- 
length is comparable or smaller than the interparticle spacing. This simple 
cut-off regularization scheme cannot be used to calculate the actual finite 
energy E of the depletion because the result is clearly cut-off dependent. 
One must either extract the low energy physics from the formally infinite 
expression by making use of effective field theory and dimensional regular- 
ization ^Ij or use the finite expression for E derived in |121ll3j via a modified 
pseudopotential. However, to evaluate the Casimir energy, it is perfectly fine 
to use (jSJ and in fact this is what was done in [7j. The reason is that the 
Casimir energy is the difference between two energies - the discrete and the 
continuum - and this difference turns out to be independent of the cut-off 
A. The Casimir energy in the BEC picks out the long wavelength behavior 
near k = and is therefore oblivious to the ultraviolet cut-off. We will show 
explicitly that the terms in the Casimir energy correspond to terms in the 
series expansion of E{k) about k = 0. 

We now evaluate the Casimir energy density E casimir- This is equal to the 
difference between the energy density of the discrete and continuum modes 
(both bulk and surface terms). As in 0, we evaluate the Casimir energy 
via the Euler-Maclaurin formula |15j : 



Ecasimir = V /(n) - / /(n) + - /(O) 



n=\ 

oo 



S i|?T/*-'(0) = - |y(0) - § /'(O) - f /=(0) + . . . 

p=l ^ ' 



(6) 



90 d^ +315(2/z)V2d5+^^^ ^ 

90 d^ 315 vd^ ^' ' 

where v = is the speed of sound in the BEC and /^^~^(0) are 

odd derivatives evaluated at zero. The leading order result, Eg^aiar = 
—TT^ v/{90d^) = — 0.10966u/d^, is exactly equal to the Casimir energy den- 
sity of a massless scalar field obeying a linear dispersion relation and confined 
to parallel plates with periodic boundary conditions ^1 EI • The massless 
scalar field propagates with wave velocity v equal to the speed of sound in 
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the BEC instead of the speed of hght. The next term is a 'Bogohubov' cor- 
rection that arises because the Bogoliubov dispersion relation is nonhnear 
and contains higher powers of k when expanded about k = 0. The magni- 
tude of the ratio of the correction Esogo = 2 7r^/(315 w d^) to the leading 
order result Egcalar is much smaller than unity since the plate separation d 
is assumed to be much greater than the healing length i.e. 



scalar] A'' d 

Therefore, to a very good approximation, the Casimir energy of the BEC 
corresponds to quantum fluctuations of a massless scalar field propagating at 
the speed of sound (i.e. acoustic phonons). This correspondence is not only 
significant conceptually but also computationally. Results on the Casimir 
energy of massless scalar fields can be applied to the BEC to obtain leading 
order terms. We apply this principle in the next section to obtain finite-size 
corrections to the parallel plate scenario. 

It is worth noting that the terms in the Casimir energy correspond to terms 
in the series expansion of E{k) about k = i.e. 

V"" . .0 k^ 



E{k) = k ^k^ + 2iJL-k^-ii = -— + vk-k^ + — -—^k^ + 0{k'^) . (8) 

In other words, the Casimir energy picks out the long wavelength behav- 
ior near k = term by term. The first term in the above expansion is a 
constant and does not contribute to the Casimir energy. The next term pro- 
portional to k is the linear dispersion relation of a massless scalar field and is 
responsible for the leading order Casimir energy Egcalar- The next term, pro- 
portional to fc^, does not contribute to the Casimir energy. This can be seen 
from the fact that k'^ = (2n7r/d)^ + contains a term proportional to 
whose odd-derivatives evaluated at zero are zero or alternatively, from zeta 
function regularization we obtain C(~2) = . The k"^ term is responsible 
for the correction Esogo- Let us show this explicitly. The function /(n) for 
/fc3 = (4n27rVd2 + r2)'/' is /(n) = 1/(2 vr) (4n2 Tr^/d^ + r^)'/' r dr. The 
odd-derivative terms /^^^^(0)/(2p)! are all zero except — 56/^(0)/6! = 
47r''/(315d^). Multiplying by the A;^ coefficient 1/(2 yields 27r''/(315i; d^) 
which is equal to Esogo- 

The Casimir pressure is readily obtained by taking the derivative of the 
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Casimir energy density 



1. 



casimir 



^-^casimir \ 7 TT V TT 



+ TT^^ + O (^"' d"') • (9) 



ad Jn 180 d^ 35vd 
The leading order Casimir pressure is 

^scalar — -j^gg ^4 • l-^u; 

It is negative (attractive) and inversely proportional to the fourth power 
of the distance as in Casimir's original calculation for the electromagnetic 
field p. Casimir obtained — vr^ c/(240 d^) for the attractive force per unit 
area between parallel conductors in vacuum. The Casimir pressure in the 
BEC is considerably weaker than in the electromagnetic case because the 
speed of sound v is orders of magnitude smaller than the speed of light 
c. Although Casimir's original calculation was performed in 1948, one had 
to wait nearly 50 years before the Casimir force between two conductors 
was conclusively confirmed by experiments starting with Lamoreaux's 1997 
landmark experiment with a torsion pendulum (2| and then by Mohideen 
et al. 1998 experiment with an atomic force microscope Here also, 
one can expect that theory is considerably ahead of experiments and that 
measuring the much smaller Casimir force in a BEC is something for the 
next generation of BEC experiments. 



3 Leading order Casimir pressure in three dimen- 
sional cavity and finite-size corrections 

In this section we obtain expressions for the leading order Casimir pressure in 
a zero-temperature homogeneous dilute weakly-interacting BEC confined to 
a three-dimensional rectangular cavity of arbitrary dimensions Li x L2 x L3 
with periodic boundary conditions (all three lengths are assumed to be much 
greater than the healing length). We then determine finite-size corrections 
to the parallel plate result Egcaiar and Pscaiar obtained in the last section. 
As in the previous section, the only parameter of the BEC that enters in 
our calculations is the speed of sound v = (2^)^/^. We take advantage of 

^Note that the derivative of the speed of sound v with respect to d (by definition keeping 
the number of atoms A*' constant) is not zero but equal to —v/{2d). If the derivative had 
been assumed to be zero, the numerical factor of the leading term in @ would have been 
be -6/180 instead of -7/180. 
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the equality between the Casimir energy of a massless scalar field and the 
leading order Casimir energy of a BEC by making use of recently derived for- 
mulas for the Casimir energy of massless scalar fields propagating with speed 
u in a d-dimensional rectangular cavity of arbitrary lengths Li, L2, . . . , 
[T7j . The speed v in this context is the speed of sound in the BEC. The 
Casimir energy in was conveniently expressed as a compact analytical 
part plus remainder. For periodic boundary conditions (inserting h) the 
Casimir energy is given by |17j : 

Rj is a remainder term expressed as an exponentially fast converging sum: 
00 <x> 2(nL,+i)^ K (|^^(£,Li)2 + ... + (£,.L,)2 ) 

R, = x; E — —1 • 

„=1 «,=-oo TT [{ii Li)2 + + LjY\~ 

(12) 

The prime means that the case where all ts are zero is excluded. There 
is no remainder for j = (it starts at j = 1). As discussed in |17j . the 
remainder term is small if we label the longest length Li, the next greatest 
length L2, etc. In three dimensions, we therefore label the lengths such that 
Li > L2 > L^. There is a clear physical interpretation to the analytical 
and remainder part in which is discussed in section 3 of ^7j- The 
analytical part is the sum of individual parallel plate energies out of which 
the rectangle is constructed while the remainder is a small contribution 
due to the nonlinearity of the energy. For a cube with periodic boundary 
conditions the remainder is 1.5% of the Casimir energy (a table of numerical 
results can be found in section 5 of JZl)- For three arbitrary lengths, the 
remainder is even smaller than in the case of the cube. 

From (|llj) we easily obtain the Casimir energy in three dimensions {d = 3) 
for periodic boundary conditions: 

(13) 

The remainders Ri and R2 are sums over modified Bessel functions given 
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by JUl) i.e. 



n=l e=i 



An L2 



Li 
L2 



R2{Li,L2,Ls)=Y^ 



00 00 2n^/^ 7^3/2! 2 7rn 



L 



is 



rt=l £1/2=— 00 



vr 



+ 



3/4 



(14) 



Note that only the ratios of lengths appear in the remainders and that 
Li/L-^ > 1, L2/L3 > 1 and L1/L2 > 1 since the lengths are labeled such 
that Li > L2 > L-^. The Bessel function Ki{2 it n I Li / L2) has a maximum 
value of 9.87 x 10~^ which occurs when n = i=l and L1/L2 = 1. As the 
ratio of lengths increases, the Bessel function decreases exponentially fast 
and becomes tiny quickly. For example, if L1/L2 = 10, then Ki[2tt 10) is of 
order 10~^^. One obtains the same order of magnitude for 7^3/2(2 vr 10). 

In the parallel plate scenario, the plate separation is the smallest of the three 
lengths and is therefore -L3. The Casimir energy per unit area is then 



1 



L1L2 



hv^ _ 7r2 _ C(3) /L3 

90 27T \L2 



Q\Li) L2 



where R is the remainder contribution given by 



R 



-vr 



i?i(Li,L2) - 7ri?2(Li,L2,L3 



(15) 



(16) 



We recognize the leading term as the parallel plate result 
^scalar = "'''"^ ^7;/(90d^) with d = L3. The other three terms are the finite 
size corrections to the Casimir energy. Let us now calculate the pressure 
along the plate separation L3. This is given by the negative derivative of 
the Casimir density with respect to L3. Keep in mind that the velocity v has 
a dependence on L3 and that dv/dL^ = —v/{2L^). The Casimir pressure 
on the plates is given by : 



P 



1 dEp 
-Li Lo d L-i 



111 
180 



hv / 

It V 



(17) 
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Pr is the contribution of the remainder given by 

PR = -^{^yRi-'^-^R2 + 7TL,R', (18) 

where the prime is a derivative with respect to L3. These derivatives are 
trivial to evaluate using the expressions H14|) and yield again Bessel functions. 

The interpretation of H17|) is straightforward. When Li ^ 00 and L2 — > 00, 
only the first term survives and one recovers the leading order parallel plate 
result Pscaiar = — 7 TT^ ^ z;/(180 d^) with d = L3. The other three terms in 
H17|) are the finite size corrections and depend on the ratios L3/L1, L3/L2 
and their inverse. Pr is orders of magnitude smaller than the other two 
correction terms because, as we have seen, the Bessel functions and their 
derivatives are tiny even in the case where the lengths Li and L2 are equal to 
L3. It is therefore an excellent first approximation to drop Pr to determine 
the finite size corrections. Let us therefore make a few quick numerical 
calculations. We quote results in units oi hv/L'^. In these units Pscaiar = 
-77rVl80 = -0.383818. Let us begin with the extreme case where Li and 
L2 are equal to L3 i.e. a "cube". Then P = -77rVl80-C(3)/(47r) -7r/12 = 
—0.741274. The force remains attractive but is much stronger. It constitutes 
a 93% difference from Pscaiar- Clearly, forces depend on the finite sizes of Li 
and L2. We are however interested in finite-size corrections to the parallel 
plate scenario where Li and L2 are at least a few times longer than L3. 
Consider then the case where Li and L2 are 5 times longer than L3. Then 
P = —0.386678. This constitutes only a 0.75% change from Pscaiar- If -^1 
and L2 are 10 times longer than L3, the pressure reduces to P = —0.384175, 
which constitutes only a 0.09% change from Pscaiar- Therefore, for any 
realistic parallel plate scenario where Li and L2 are at least a few times 
longer than L3, the Casimir pressure is dominated by Pscaiar - 

The result for the pressure given by ()17() . including the remainders, is a 
leading order result. It does not include the next order, the Bogoliubov 
corrections which arise from the term in the series expansion of E{k) 
given by We saw last section that these Bogoliubov corrections are 
small. Nonetheless, if we wish to include them, one needs to obtain the 
analog formulas to (|llj) for higher power dispersion relations and this is 
work for the future. However, we can already include the dominant Bo- 
goliubov correction. There is a separate Bogoliubov correction associated 
with each of the four terms in (|17() . If we label the four terms Pi,P2;-P3 
and P4 then there exists a Bogoliubov correction for each term which we la- 
bel Pin ,P2h iPsh and P4„ . The dominant Bogoliubov correction. 
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^iBogo ' is associated with the leading term Pi = —7iT'^hv/{180L'^). Its value 
has already been calculated in @ i.e. Pig^^^ = 7r^/(35f d^). To match the 
units and notation of this section, we insert fi and m and let d = L3. This 
yields Pig^g^ = ti"^ tt^ / (140 m? v L^). So we can already include the most 
important Bogoliubov correction into the result H17|) by replacing P with 
P + Pig^g^ ■ Pisogo is clearly not a significant correction to Pi , but it can 
still compete with the finite-size corrections P2 and P3 when Li and L2 
are sufficiently larger than L3. Of course, both corrections are then small. 
To summarize, in the parallel plate scenario Pi = Pgcaiar is the dominant 
Casimir pressure in the BEG and corresponds to the parallel plate result 
for a massless scalar field. There are then small finite-size corrections to 
Pscaiar- If One is interested in keeping track of tiny corrections, Bogoliubov 
corrections become comparable to finite-size corrections when Li and L2 are 
sufficiently larger than L3. 
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